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Abstract Ionizing radiation is known to delay the cell cycle progression. In particular 
after particle exposure significant delays have been observed and it has been shown 
that the extent of delay affects the expression of damage such as chromosome aberra- 
tions. Thus, to predict how cells respond to ionizing radiation and to derive reliable 
estimates of radiation risks, information about radiation-induced cell cycle perturba- 
tions is required. In the present study we describe and apply a method for retrieval of 
information about the time-course of all cell cycle phases from experimental data on 
the mitotic index only. We study the progression of mammalian cells through the cell 
cycle after exposure. The analysis reveals a prolonged block of damaged cells in the 
G2 phase. Furthermore, by performing an error analysis on simulated data valuable 
information for the design of experimental studies has been obtained. The analysis 
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showed that the number of cells analyzed in an experimental sample should be at least 
100 to obtain a relative error less than 20%. 

Keywords ionizing radiation • cell cycle delay ■ Monte Carlo simulation 



1 Introduction 



It is well-known that cells respond to DNA damage by activating checkpoints that delay 
the cell cycle transition in particular from Gl to S phase an d from G2 t o M p hase. 
These delays are assumed to allow additional time for repair I Li and Zoul 1 2005h ). 

After sparsely ionizing radiation such as X-rays or 7-rays mild perturbations of the 
cell kineti cs have been observ ed corresponding to the delay of about 1 hour per 1 Gy of 
exposure (jPurrot et a3 (|l980h ). In contrast, after particle irradi ation dramat i c cell cycle 
delays have been measured lasting up to 3 cell generation times dScholz et all ( 19941 ) and 
references therein). Moreover, a rapid desynchronizat ion of initial l y syn c hronous cell 
populations has been observed after particle exposure ( Scholz et all ( 19941 ) ; iRitter et all 
(|l99dl2000h . The differences between the effects induced by both radiati on qualities 
may b e directly attributed to spatial differences in the energy deposition ( Ritter et all 
( 200d ); lGudowska-Nowak et all i|2005l )). 

Cell cycle delays as observed after radiation exposure are important for the in- 
terpretation of other biological experiments. For example, recent reports have shown 
that cell cycle delays interfere with the time-course of aberrations visible in cells at 
the first mitosis post-irradiation. In particular, after particle exposure which delays the 
cell cycle progression more than sparsely ionizing radia tion, a dras t ic increase in the 
aberration yield with sampling time has been observed IjRitter et a 3 (|l996l . l200cl )) and 
it has been recently shown that the average time to enter the first mitosis correlates 
directly with the aberration burden of a cell l|Gudowska-Nowak et a 1 (2005)). In other 
words, cells entering mitosis at later times harbor more aberrations than those entering 
mitosis earlier. Since the frequency of aberrat ions expressed in first cycle me taphases 
is used to determine the absorb ed dose (e.g. Cucinotta and Durante! ( 20061 )) and to 



derive cancer risk estimates (e.g. iMateuca et all 1 20061 ) ). the damage sustained by the 
whole cell population has to be determined. This can be archived by the analysis 
of samples collected at multiple sampling times covering the whole interval from the 
first to the last cells reaching mitosis. Then, the total yield of ab errations can b e deter - 
mined by a mathem atical approach (i.e. integration analysis, see Kaufman et a] (|l974l ): 
IScholz et all (|l998l )). 

Experimental and theoretical studies of the cell cycle progression a r e also of primary 
i ntere st in the cont e xt of c a ncer therapy | Hahnfel dt and Hlatkvl dl996l ): lM"ontalenti et all 
( 19981 ); lErba et all J2002I ); iBasse and Ubeziol (|2007l ); I Wilson! (|2007l )). since the thera- 
peutic response of solid tumors is known to depend not only on the administrated dose 
of ionizing radiation or chemotherapy agents but also on repopulation and redistribu- 
tion of cancerous cells. Cancer therapies may target specific phases of the cell cycle 
by blo cking or delaying the progress through one or more phases (see lMontalenti et all 



The specific aim of our project was to further elucidate the complexity of particle- 
induced cell cycle arrest. Investigations into the effects of particles are becoming in- 
creasingly pertinent in light of rapidly growing interest in this type of radiotherapy 
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( Amaldi and Kraft! (|2007h ). Furthermore, for the planning of manned missions to Moon 
and Mars, a better knowledge of the action of charged particles is needed, since the 
main contribution of dose during a space mission outside the magnetic shielding of the 
Earth originates from galactic cos mic rays, which are heavy part icles from the most 
frequent protons to up to iron ions ( Cucinotta and Durante! 1 20061 )). 

In this study, we present a method for retrieval of information about all previous 
phases of the cell cycle, when only the mitotic index of the experimental cell sample 
is known. We use the met hod to analyze the data measured for mammalian V79 cells 
after exposure to Ar ions i Ritter et all ( 2000l )) . V79 cells re present a frequent l y use d 
model system to study genotoxic effect s of ion i zing radiation (IRitter et a 1 (1996, 2000); 
IWevrather et all ill999l ); lGroesser et all l|2007i ); |Pathak et all < |2007h ) or chemical agents 
( Virgilio et all ( 20041 ) ). Information on the cell cycle progression was available in the 
form of subsequent measurements of the mitotic index in control and irradiated s amples 
(Fig. 2) and duration times of cell cycle phases of control cells ( Sinclair! 1 19691 ); IScholj 
(2003)). Since in the experiment not only the mitotic index, but also the aberration 
yield has been measured, the approach may allow us to study in future the progression 
of cells carrying a different number of aberrations. 

Mitotic index is the measure of the number of cells undergoing mitosis at a given 
time. As the duration of cell cycle phases varies from cell to cell, the mitotic index 
depends on the duration distribution of mitosis, but also on the corresponding du- 
ration distributions of the previous phases. We retrieve these distributions using a 
multi-dimensional fit based on general ass umptions concern i ng th e distribution shape, 



which are inferred from other experiments ( Montalenti et all . Il998l ). In this way, having 
only the mitotic indices for control and irradiated cells, we were able to compare their 
progression through the whole cell cycle and predict which phases are most vulnerable 
to irradiation. The fitted parameters of the duration distributions of cell cycle phases 
have also been used to estimate the experimental error of the mitotic index measure- 
ment, depending on the number of cells in the sample. The estimation has been carried 
out using the Monte Carlo simulation of the cell cycle progression for a given number 
of cells. 



2 Materials and methods 

2.1 Experimental data used for the analysis 

To examine in more detail charged particle-induced cell cycle progressi on delays in 
the fir st post-irradiation cycle, previously published data were reanalysed ( Ritter et all 

i1qq3)). 



For the experiment, V79 Chinese hamster cells were synchronized by mitotic shake 
off, a method based on the selective detachment of mitotic cells from monolayers by 
shaking. Mitotic cells were plated in Petri dishes and about 2h after seeding, when the 
cells had attached and divided and progressed into Gl-phase, the exposure was done. 

Irradiation with 10.4 MeV u _1 Ar ions (LET = 1226 keV/im -1 ) was performed with 
a fluence of 10 6 particles cm -2 corresponding to a dose of 1.96 Gy, respectively. At the 
time of irradiation, at least 95% of cells were in Gl phase as determined by flow cytom- 
etry. Immediately after exposure, 5'-Bromo-2'-deoxyuridine (10 /img ml -1 ) was added 
to the samples to distinguish between metaphases in the different post-irradiation cy- 
cles. Cells were harvested at multiple sampling times covering the time interval of 2 to 
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Fig. 1 Schematic representation of the cell cycle consisting of 4 distinct phases, namely Gl, 
S, G2 and M. For V79 cells used in the present study the mean dur a tion times of the t hese 
phases are 2.25 h, 6.5 h, 1.5 h and 0.75h, respectively (ISincIairl dl969T ); IScholz et all [U994T 0. 



3 cell generation times. Each sampling time was preceded by a 2h colcemid treatment 
(0.1 A*mg ml -1 ) to accumulate mitoses. Since the main purpose of the experiments was 
the analysis of the effect of cell cycle delays on the expression of cytogentic damage, 
chromosome preparations were made according to the standard procedure and slides 
were stained with the fluorescence-plus-Giemsa technique. At each sampling time chro- 
mosomal damages was scored in 100 first cycle metaphases IjRitter et a ] fcpOCh ). To 



gain information on the cell cycle progression, the mitotic index was determined for 
each sample by the direct scoring of 1000-2000 cells on the slides and the cell gener- 
ation of at least 200 metaphases was recorded. Due to experimetal limitation (access 
to the particle beam and number of cells that can be synchronized), more detailed 
measurements of the cell cycle progression, with a better time resolution, for example 
by flow cyt ometry, were not f easible. A detailed description of the experimental setup 
is given in 



( Ritter et all (|2000h ). 



The passage through the cell cycle has been simulated by means of a kinetic model 
in which the duration of each phase is taken as a random variable characterized by 
its mean and dispersion (see below). The motivation for such an approach is the rapid 
radiation-induced desynchronization of initially synchronous cell populations. The cell- 
cycle time (the time interval between cell divisions) becomes then a random variable 
whose statistical properties can be inferred from the analysis of the frequencies of cells 
observed in different phases at different times. However, as mentioned above, the avail- 
able information on the cell cycle kinetics relies solely on the measured mitotic indices. 
Therefore, for the analysis, these values were used to deduce (a posteriori) duration 
times in all phases before the mitosis. In the following sections, we aim to expand this 
approach and to present a model which is general enough to be applied not only for 
the analysis of radiation-induced cell cycle delays but also for the an alysis of multi- 
compa rtment cell populations perturbed by other cancer therapies ( Kohandel et all 



(2007)). 
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2.2 Mathematical model of cell cycle progression 

In many experiments a positive correlation between radiation dose and the dura- 
tion of cell-cycle delays was found. Although such findings were usually quantified in 
terms of a linear relationship betwe en phase duration and dose i Zaider and Minerbol 



( 1993h : lHahnfeldt and Hlatkvl (jl996h ). a more detailed analysis points to a direct cor 



relation between cell cycle de lay and the number of aberrations carried by a cell 
( Gudowska-Nowak et af (|2005l) ). This effect is responsible for the loss of synchrony 



of the population and can be illustrated by interpreting the (normalized) mitotic index 
as a frequency histogram of times spent by cells before the actual division happens. The 
cell cycle kinetics can be then investigated by treating the durations of the four phases 
as independent stochastic variables ha ving probability den s ity fu n ctions described 
ter ms of two adjustable parameters ( Zaider and Minerbol ( 19931 ); iMontalenti et all 



(Il998h . The parameters of such distributions are fixed in time (stationarity of the 
distribution is assumed), whereas the choice of a particular probability distributi on 
function has been noted not to be critical for the final result ( Hartmann et all ( 1975h ). 



For V79 Chinese hamster cells, as used in our simulations, the mean duration 
times of the c e ll cyc l e phases are tm = 2.25 h, tg — 6.5 h, t<j2 = 1.5 h, t^i = 
0.75h l|Sinclairl l| 19691 ): I Scholz et 3 (|l994l ^. The phase duration for an individual cell 



is given by a certain probability distribution _Dp/j(r), where r is the time which a 
given cell had already spent in the current phase. Consequently, a single cell of phase 
age ph is assumed to leave its current phase within a time interval [r, r + dr] with a 
probability 7?p/ l (r)cir. The mean number of cells dN p h^{t) leaving the phase ph within 

an infinitesimal time interval [t, t + dt] is defined by the mean flux dN p^W 

The experimental input to our model is given by the mitot ic indices scored f or 
exposed and control cell populations at subsequent 2h intervals IjRitter et a ] (|2000l)). 



Experimentally measured values of the mitotic index MI(i) at a timestep i 

are defined as the increase in the total number of cells in relation to their number at 
the previous sampling time. Since the time until all cells reached the first mitosis was 
quite long compared to the average cycle length (Fig. [2j , c hanges in the p o pulation size 



quite long compared to the average cycle length (t lg. [zj , changes m the population size 
due to cell division have to be taken into account ( Kaufman et all ( 19741 ); IScholz et all 
(HI)). 



To determine the increase in the number of cells in relation to their initial number 
one expresses then the corrected mitotic index in the form 

/•x ,^ T /.s^(i-l) N(i) — N(i — 1) 
MI corrected (t) = MI(i) i— > = 1 ; I '- (2) 

In turn, the summed mitotic indices up to the time step n 

S{n) = J2 MI corrected (i) = N(n) N ~ N ° (3) 

give then the fraction of cells which have already completed a full cycle, i.e. went 
through mitosis up to the given time. The experimental results show that this quantity 
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Fig. 2 Corrected mitotic in dex of V79 cells expose d in G i phase to 1,96 Gy Ar ions (10.4 
MeV/u): experimental data jGudowska-Nowak ct al (2005); Rittcr ct al (2000)) along with 
fits and simulations of No = 1000 cells are shown. The mitotic index curve was fitted with 
a lognormal probability density function with parameters fi and a (cf. Eqs.Jg] 0) and the 
quality of the fit was determined by using the standard x 2 test which yielded almost the same 
results in two cases. However, the underlying cell cycle kinetics is quite different for Fit 1 and 
Fit 2 (see Sec. [3] and Figs. [4] [5J. In Fit 1, the Gl phase is largely dispersed whereas other 
phases have well defined duration times of negligible variance. In Fit 2 the phases S and G2 
are dispersed and Gl and M duration times are assigned constant values. 



tends to a constant level, which yields approximately 100% for control cells and less 
than 100% for irradiated cells, as displayed in Fig. [3] 
The mean 

fl ux gjWi?) defineg 

a mean number of cells dN p ^_, (t) leaving the phase 
ph within an infinitesimal time interval [t,t + dt]. The term "mean" is understood here 
as an average over a large number of identical experiments, each with an initial number 
of iVo cells. Note that t stands here for the "absolute" time, different from the phase 
age r. In order to relate the description in terms of a mean 

flux dN vh ^(t) baged 

on 

the ensemble approach to the description of the time evolution of a single cell whose 
probability of leaving the current phase is Dp^(r), we assume that we deal with a large 
ensemble of samples, each initially containing iVo cells in the same phase ph. Moreover, 
the sample populations are assumed to be of the same phase age, e.g. r = 0. We start 
our "numerical experiment" at time t = 0, so that t = t, and count the cells which 
leave the phase within the age r and r + dr. Their mean number will be given by: 



dJV ph-» = N F ph (r)dT, (4) 

where F p ^ (r) is a probability distribution which gives the likelihood of a cell leav- 
ing the phase at the time t = r. 

On the other hand, if we trace the behavior of a particular cell, we can define D p h (t) 
as the conditional probability distribution of leaving the phase at age r, provided that 
the cell had not completed the phase earlier: 
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Fig. 3 Cumulative fraction of cells which left the M phase of the first cycle, compared with 
the experimental results. For the explanation of fits, see Sec. [3] 



Fphjr) 
JoF ph (r')dr'- 



D ph(r) = r^TTTmr- (5) 



Function Fpi^t) in the above formula represents the distribution of phase dura- 
tion times in an initially synchronized population of cells. We follow the suggestions 
of Zaider and Minerbo (1993) and Montalenti et al. (1 998) and assume a lognormal 
distribution of phase duration ( Engen and Landel ( 1996h ): 



tJ2naL \ 2a P h 



pit 

with parameters and a p ^ determining the mean 



and variance 



E[t] = e^ h+ -^ (7) 



E[(t - E[t]f] = e 2 ^ h+a p h (e a ^H _ X ). (8) 



In the above expressions, fi and a are the parameters to be fitted. 

It should be noticed that a log-normal distribution is not the only possible choice 
to be postulated as suitable for the description of the random distribution of phase 
duration times. The other possibility might be, e.g., the gamma distribution res ulting 



from a mixture of exponential distr i butions expected in simple ren ewal processes (Feller 



Irom a mixture ot exponential distributions expected m simple rene^ 
( 19681 ); IZaider and Minerbol (|l993l l; lHahnfeldt and Hlatkvl (|l996T l) 



Ensemble analysis of the mean flux allows splitting dN^^it) into a sum of con- 
tributions from subsequent cycles: 

dN ph ^(t)^J2 dN pL( t ) (9) 
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dNgJ(t) is defined as the product of the probability of leaving the Gl phase within 
the time interval (t, t + dt) and the mean total number of cells at time t: 

dN^(t)=dP G1 ^(t)N(t), (10) 

which means that 

dP G1 ^(t) = F Gl {t)dt. (11) 

In the first cycle N(t) = Nq, so 

dN^(t) = dP$_>(t)N = N dtF G1 (t). (12) 

The probability of leaving the S phase in the first cycle at a certain time will depend 
on the probability of leaving the previous phase, and thus 

dN i s 1 \ f (t) = dP { s 1 l(t)N a = N a dt F G1 (t-T)F s {r)dT = N dtFi(t), (13) 

Jo 

and, analogously: 

<W%L(t) = dP { G l^{t)N = N dt f* F x (t - r)F G2 {r)dr = N dtF 2 (t), (14) 

JO 

dN^(t) = dP^(t)N = N dt [ F 2 (t - r)F M (r)dT = N dtF 3 (t). (15) 

Jo 

In the moment of leaving the M phase, cells divide: two Gl cells are produced and 
therefore, the number of cells in Gl phase (and thus the number of those leaving Gl 
phase) is twice larger: 

dN G 2 }^(t) = 2dP&(t)N = 2N a dt [ F 3 (t - T)F G1 (r)dT = 2N dtF 4 (t). (16) 

JO 

In this manner we generate consecutive contributions to Eq. ([9]). Knowing the mean 
flux as given by Eq. ((9} we can calculate other mean quantities of interest, such as 
mean number of cells found in the phase ph at a given time after irradiation. This 
method gives therefore a global view on the cell cycle kinetics. 

In order to reconstruct the fate of an individual cell, we have used the Monte Carlo 
method. The computer simulation tracks the cycle of each single cell separately. Ini- 
tially we generate an ensemble of TVo cells. Each cell is assumed to start in the Gl 
phase and its phase age is r = 0. In subsequent time steps each cell goes through its 
individual cycle, and in every time step dr each cell has a probability _D p / l (r)dr to 
leave its current phase. T he algorithm is built based on the following scheme: 

1. The cell is in the phase ph. Its phase age is r. 

2. A random number x, < x < 1 is generated. 

3. If x < D p ^(r)dT, the programme proceeds to 4., else it goes to 6. 

4. If the next phase "ph + 1" 7^ Gl, then the cell exits from its current phase: 
ph = "ph + 1", t = 0. Go to 1. 

5. If the next phase "ph + 1" = Gl, then the cell exits from its current phase and 
replicates: ph — Gl, t = 0. Additionally, a new cell is generated, whose current phase 
is Gl and r = 0. The cycle starts with going to 1. 

6. The cell remains in its current phase: r = r + dr. At the next step the programme 
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Fig. 4 Control cells: Probability density functions for leaving a phase in phase age r (top 
panel: Fit 1, bottom panel: Fit 2) 



starts with 1. 



All integrations to compute D p h (r) from ([5]) have been performed in our simulation 
using the Euler method with discrete time steps At. 



The above described approach allows to examine the kinetics of a single experiment. 
Computing the variance of quantities of interest, we can estimate how much a single 
realization of the experiment deviates from the ensemble average. 
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Fig. 5 Irradiated cells: Probability density functions for leaving a phase in phase age r (top 
panel: Fit 1, bottom panel: Fit 2), cf. Eqs. lT8l[T9t . 



3 Results 



Using the methods described above, we have first fitted the experimental results, i.e. 
the mitotic indices (see Sec. l2.ip . Since the experimentally accessible information refers 
to cells initially synchronized in Gl phase, we assumed an "ideal synchrony" for the 
numerical analysis , i.e. at time t = all of cells were assumed to be of the phase age 
r = 0. 
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Control cells 




Fit I 


Fit II 


Phase 


Median 


Mode 


CV 


Median 


Mode 


CV 


Gl 


1.806 


1.164 


1.346 


2.25 


2.25 




S 


6.500 


6.500 


500 


6.467 


6.401 


9.876 


G2 


1.5 


1.5 




1.000 


0.445 


0.895 


M 


0.75 


0.75 




0.75 


0.75 




Irradiated cells 




Fit I 


Fit II 


Phase 


Median 


Mode 


CV 


Median 


Mode 


CV 


Gl 


2.250 


2.250 


250 


2.398 


2.398 


250 


S 


6.500 


6.500 


500 


6.467 


6.401 


9.876 


G2 


2.744 


0.146 


0.237 


2.389 


0.084 


0.192 


M 


0.75 


0.75 




0.75 


0.75 





Table 1 Median, modal values and coefficients of variation for the fitted phase duration time 
distributions Fp^. 



3.1 Control cells 

The first step in our analysis was fitting the parameters of the F p h distributions (0 
for each phase, in such a way that they together made up the time-course of mitosis 
(|15[) which fitted the experimental corrected mitotic indices ([2]). 

We performed the integrations numerically with a time-step At = 0.03 h (Eulcr 
method), obtaining an approximate flux: 

AN ph ^(t) _ dN ph ^(t) 

At ~ dt [ ' 

Having evaluated the outflow from the M phase, we could further easily compute the 
total mean number of cells at a given time, from which the (mean) corrected mitotic 
index ([2]) was constructed. (We recall that the word "mean" is understood here as an 
average over a large number of identical experiments.) 

Since we would have had to fit up to 8 parameters, we decided to simplify the task. 
The results cited in Sec . 1 2 . 1 1 delivered the values of mean duration times of each phase: 
Jqi = 2.25 h, tg = 6.5 h, tQ2 ~ 1.5 h, ijtf = 0.75h. Therefore, we used them as the 
four (frozen) parameters needed for modelling and focused on four others (we chose 
(Tph, which can be treated as a measure of fitted distribution's width) which had to 
be chosen by the best fit. The initial number of cells assumed for the modelling was 
iVo = 1000 and the parameters were fitted by using a Metropolis algorithm. 

We found that the result of the fitting (see Fig. [2]| was ambiguous and there were 
two possible least-square fits with \ 2 value of the same order. The inferred coefficients 
displayed a significant variability between both data sets (see Fig. [4] and Tab.[T]): 

Fit 1: 

a G1 = 0.663, a s = 0.002, a G2 = 0.000, a u = 0.000, \ = 55.0 (18) 

Fit 2: 



a G1 = 0.000, a s = 0.101, a G2 = 0.900, a M = 0.000, x = 50.0 (19) 
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Fig. 6 Flux from phases Gl and S (number of cells passing per lh) at time t, compared for 
control and irradiated cells (Fit 2). 



Although the two possible cell cycle schemes produced almost identical mitotic 
index curves, the intra-cycle distributions differed very much. In Fit 1, the Gl phase 
is largely dispersed whereas other phases have well-defined (deterministic) duration 
times. In contrast, in Fit 2, S and G2 phases are dispersed with the duration times of 
Gl and M phases being well specified. 



3.2 Irradiated cells 



In order to reproduce the mean corrected mitotic indices of the irradiated cell popula- 
tion we performed a similar fitting procedure as mentioned above. To proceed, we made 
the biologically justified assumption that irradiation affects only the characteristics of 
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Fig. 7 Flux from phases G2 and M (number of cells passing per lh) at time t, compared for 
control and irradiated cells (Fit 2). 



Gl and G2 phases, since checkpoints in Gl and G2 are known to block the cell cycle 
progression to give the cell time to repair and continue cycling or to undergo apoptosis 
(|Li and Zoul (|2005r i). We also expected that the mean duration times of both phases, 
as well as the variances of the duration times should increase because the damaged 
cells would be blocked in either one of them for a longer time. Therefore, we fitted only 
the parameters ten, cr Gl>^G2i °"G2> taking others as fixed and referring to their values 
obtained in the former fit for control cells. While doing so, we did not introduce any 
additional mortality parameter. Instead, we assumed solely that the damaged cells will 
stay in Gl or G2 phase for a very long time. 



For the parameters corresponding to Fit 1 the analysis yielded: 
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t Gl = 2.25, a G1 = 0.004, t G2 = 11.90, a G2 = 1.713, X 2 = 5.6, (20) 
whereas for Fit 2 it resulted in: 

t Gl = 2 AO, o G1 = 0.004, t G2 = 12.70, a G2 = 1.828, X 2 = 5.8 (21) 

In both fits the F G2 (r) distribution function became wider and its mean duration 
time strongly increased (Fig. [5] see also Tab.[T}. 

In turn, the F G i(r) distribution in Fit 1 became narrower in comparison to the 
corresponding distribution for the control cells, which implied that Fit 1 was unrealistic. 
This finding is consistent with previous flow cytometric studies showing that V79 cells 
suffer short cell cycle delays in Gl phase, but more pronounced delays in S- or G2-phase 
(Scholz et al. 1994). The low activity of the Gl checkpoint in V79 cells might result 
from mutations in the p53 gene (Chaung et al, 1997). We have therefore concluded 
that Fit 2 was confirmed as better matching the biological scenario: here the mean 
duration time of Gl increased and the distribution function became wider. The very 
long tail of F G2 (r), as predicted in this case, suggests that damaged cells may remain 
blocked in the G2 phase for an extremely long time with a non-zero probability as it 
has been experimentally shown by Scholz et al. (1994). 

Figs. [6] and [7] compare the progression of cell fluxes between subsequent phases for 
control and irradiated populations, respectively, while the number of cells remaining in 
a given phase at time t is shown in Figs. [8] and [9] In the first cell cycle a notably slower 
outflow of cells from the G2 phase is visible which corroborates with the number of 
cells remaining in that phase up to 16h after exposure. Cells blocked in that phase also 
diminish the fraction of those which are able to enter mitosis up to the same time-point. 
As expected, the cumulative number of cells which enter mitosis is much lower in the 
irradiated population than in the control sample and can be well correlated with the 
fraction of cells arrested in G2. 



3.3 Determination of experimental errors 

Another important finding of our work is that the applied approach allows for es- 
timation of inter-experimental differences resulting from the stochastic character of 
phase duration times. By performing a Monte Carlo simulation, described in Sec. 12.21 
we simulated a time-evolution of an ensemble of cells (here: No = 1000) using the 
set of previously fitted parameters. The number of cells at different moments of time 
was counted, and therefrom the mitotic index (JT| was derived. Using this method, the 
variance/standard deviation of the corresponding probability distribution function was 
estimated (see Fig. Ill[) . 

We performed a series of 100 cell cycle simulations with the initial number of 
cells Nq. The mean number of cells found in a given phase {N p ^(t)) and the variance 
a N (*) were calculated. Analyzing the mutual dependence of the relative error i* T p \ 
and (N p h) we found that the relative error does not depend on the initial number of 
cells Nq nor on the set of parameters of F p fr. It depends only on the current number 
of cells. A practical conclusion drawn from these results is that the number of cells 
analyzed in an experimental sample should be at least of order of 100 to obtain a 
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Fig. 8 Number of cells in phases Gl and S at time t, compared between control and irradiated 
cells (Fit 2). 



relative error less than 20%. If a better precision is required, a larger cell population 
has to be examined. 



4 Discussion 

In stu dies preceding this project (|Ritter et al (|l996ll2OO0h :l Scholz et H Il99l ) : lGudowska-Nowak ^Tal 
(2005)) the relationship between radiation-induced mitotic delay and expression of 
chromosome damage was shown. To further establish differences in progression through 
the division cycle of unirradiated (control) and exposed cells, we developed a method 
for retrieval of the information about the time-course of all cell cycle phases, when 
only the mitotic index of a cell sample is known. The method consists in a multi- 
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Fig. 9 Number of cells in phases G2 and M at time t, compared between control and irradiated 
cells (Fit 2). 



dimensional fit based on general assumptions concerning the duration distribution of 
cell cycle phases, which are inferred from other experiments (|Montalenti et 

The analysis of the experi mentally measured mitotic indices for mammalian V79 



cells after exposure to Ar ions IjRitter et al (2000)) revealed a prolonged block of dam- 



aged cells in the G2 phase. Additionally, we have shown that this approach is applicable 
directly to determine experimental errors resulting from the stochastics of phase dura- 
tion times. 

In future studies the quality of the model will be validated through the application 
to data sets generated for other cell lines and cell types such as human skin fibrob- 
lasts and lymphoc y tes whose cell-cy c le kinetics after irradiatio n is altered differently 
(|Flatt et all (|l998h ; iLee et all fcood ); iTenhumberg et all (|2007h V For example, it has 
been shown that normal human fibroblasts exposed in G0/G1 phase to ionizing radia- 
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Fig. 10 Number of cells found in M phase at time t . Black line: theoretical prediction obtained 
using the mean flux approach. Dashed line: a single simulation for Nq = 1000 cells. Error bars: 
standard deviation for the simulation. 
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Fig. 11 Relative error of the measurement vs. mean total number of cells analyzed. 



tion suffer a prolonged Gl arrest IjFlatt et aj (|l998l ); lTenhumberg et al jgOOg )), while 
human lymphocytes are predominantly arrested in G2 l|Lee et all i(2005l )1. Furthermore, 
the model may be used to analyze the progression of cells carrying different numbers 
of aberrations. 
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